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Abstract 

The manipulation of TO-phonon polaritons and the terahertz (THz) light field associated with 
them by means of an ultra-sound acoustic wave is proposed and illustrated by calculating the 
TO-phonon-mediated THz response of acoustically-pumped CuCl and T1C1 crystals. We show 
the high-contrast acoustically-induced change of the THz reflectivity alongside with multiple THz 
Bragg replicas, which are associated with the infrared-active TO-phonon resonance driven by the 
ultrasonic wave. The effect, which stems from phonon anharmonicity, refers to an operating acous- 
tic intensity 7 ac ~ 1 — 100kW/cm 2 and frequency z/ ac ^0.1 — 1 GHz, with possible applications in 
THz spectroscopy. 

PACS numbers: 43.35. +d, 71.36.+C, 78.20.Pa 
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Since the pioneering work by Kun Huang [1], the physics of infrared polaritons associated 
with transverse optical (TO) phonons has emerged as a well-established discipline. Recently, 
room-temperature polaritonics was implemented for processing and coherent control of the 
THz light field j3, S| • In addition to conventional infrared spectroscopy and continuous- 
wave Raman scattering experiments, broadband THz time-domain spectroscopy has been 
developed and successfully applied to characterize picosecond (ps) and sub-ps dynamics of 
infrared-active TO-phonons [4j and to visualize the polariton dispersion associated with 
these vibrational modes j^]. Furthermore, the THz polariton spectra allow to study unusual 
lattice dynamics, e.g., in ferroelectrics |(| and negative thermal expansion compounds j3, Sj- 
However, the spectrally-resolved control of the THz electromagnetic field associated with 
infrared polaritons still remains a very challenging task of far-infrared spectroscopy. 

In this Letter we propose an acoustic modulation of TO-phonon polaritons to drastically 
change their optical response in the THz band. Phonon anharmonicity, which can be large 
in some dielectric and semiconductor materials and particularly strong for soft TO-phonon 
modes, leads to the coupling between a coherent (pumping) acoustic wave (AW) and infrared- 
active TO-phonons. In this case one deals with an acoustically-induced Autler-Townes effect, 
which gives rise to spectral gaps in the THz polariton spectrum of AW-driven TO- 
phonons, and is akin to the acoustical Stark effect for excitons |9|. The AW-induced gaps 
A^v, which open up in the polariton spectrum and develop with increasing acoustic intensity 
7 ac , are due to the A^th-order resonant acoustic phonon transitions within the polariton 
dispersion branches. These forbidden-energy gaps strongly modify the optical response of 
TO-phonon polaritons and make possible the effective AW manipulation of the THz field. 

The spectrally-resolved AW control of the THz field propagation can also be interpreted 
in terms of Bragg diffraction of infrared polaritons by the pumping AW: For the first time we 
analyze the use of phonon anaharmonicity to create an acoustically-induced Bragg grating. 
In this case the contrast of the AW grating is dictated by the efficiency of the scattering 
channel a TO-phonon ± acoustic phonon (two acoustic phonons) <H> TO-phonon" for cubic 
(quartic) phonon anharmonicity. Thus the scattering of THz light is mediated and strongly 
enhanced by the TO-phonon resonance. This results in an anomalously short interaction 
length needed for formation of the Bragg replicas. Possible applications of governing infrared 
polaritons by using an ultrasonic acoustic field include frequency-tunable THz detectors and 
filters, Bragg switchers and frequency converters. 
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The Hamiltonian of infrared polaritons coherently driven by a cw bulk acoustic wave 
{k,w ac (/c)} is given by 

H = H + J2 [2m 4 I ac blb p + (ms/^e-^frt 6p_ k + H. c.) 

P 

+ (m 4 /a C e- W &J>&p-2k + H - c -)] > (!) 

with i7 the conventional polariton Hamiltonian of infrared- active TO-phonons [loj], 6 P the 
TO-phonon operator, cj ac = 27r^ ac = v s k, v s the sound velocity, and m 3 (7714 and 7714) the 
matrix element associated with cubic (quartic) anharmonicity. The macroscopic equations, 
which describe the control of THz polaritons by applying the acoustic field of an arbitrary 
profile, J ac = Jac(r,t), are 
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E(r,t) = -^^P(r,t), (2) 



' d 2 d 

— + 2 7tot^ + ^to + ^horhJac (r, *) + 4Q T o^3/ac /2 ( r ' *) cos(w ac t - kr) 



+ 40 TO m 4 4(r, t) cos(2cu ac t - 2kr) 



P(r,t) = ^E(r,t), (3) 



where E and P stand for the light field and TO-phonon polarization, respectively, Q^o 
is the TO-phonon frequency, is the polariton Rabi frequency, is the background 
dielectric constant in the infrared, and 7x0 is the damping rate of TO-phonons, mainly due 
to their decay into short- wavelength acoustic phonons. Equations (El)-® refer to simple 
cubic lattices with spatially isotropic long-wavelength anharmonicity and optical response. 
For cw acoustic excitations, when 7 ac (r,i) = I ac = const., Eqs. (El)-® yield the same quasi- 
energy spectrum as that of the quadratic Hamiltonian (GD). If / ac = 0, Eqs. (ED-® reduce 
to the standard TO-phonon polariton equations The forth term on the left-hand side 
(l.h.s.) of the polarization Eq. ® yields a Stark shift of the TO-phonon frequency, which 
is oc 7 ac and is associated with the quartic phonon nonlinearity. The last two terms on the 
l.h.s. of Eq. ® give rise to the Bragg spectrum of TO-phonon polaritons driven by the AW. 

The matrix elements associated with cubic and quartic anharmonicity are = 
6[V(^)] 1/2 ^ 3 (k,p-k,-p), m 4 = 12[VWfc)][V4(k,k,p-2k, -p) + V 4 (k,k, -p,p- 
2k)], and rh 4 = 24:[v /(h 2 v 2 fc)]V4(k, — k, p, — p), with v a volume of the primitive cell. The 
potential V 3 is given by V 3 (k, p - k, -p) = (l/6)[frV( 8w ac^To)] 1/2 $ (3) ( k ; P - k > ~P)> where 
$( 3 ) i s proportional to the Fourier transform of the third-order derivative of the interatomic 
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FIG. 1: (Color online) Schematic of optical excitation and backward Bragg scattering of an 
acoustically-driven THz polariton. Vertical stripes symbolize the propagating bulk acoustic wave. 

potential. There is a similar expression for the forth-order potential, V4 oc (the explicit 



formulae for V^4 and $( 3,4 ) are given, e.g., in Refs. (ill . 121]). In the long- wavelength limit 



k,p <C l/a (a is the lattice constant), relevant to the optics of TO-phonon polaritons 
$( 3 ' 4 ) are well-approximated by $^ 3 ^ = (kao)Cs and $^ = {ka^C^ [13]. The anharmonic- 



or 



ity constants C34 can either be calculated by using modern ab initio methods 
evaluated from experimental data available for some anharmonic crystals The above 

approximation of $( 3 ' 4 ) also leads to the fc-independent matrix elements m 3 and rh^ = 
in Eqs. (pQ) and ©. 

Recently, the third-order coupling constant V3 was inferred for some zincblende-type 
semiconductors (GaP, CuCl, CuBr, and /3-ZnS) which exhibit strong and dominant cubic 
anharmonicity jl6l.ll7l]. The used experimental data refers to the decay of a long- wavelength 
TO-phonon into short-wavelength longitudinal (LA) and transverse (TA) acoustic phonons. 
To adapt the inferred values of V3 to m 3 in Eq. (pQ) and Eqs. (Ej)-([3]), we use the Leibfried- 
Ludwig approximation [191 ]. The parameter m^lll 2 , which controls cubic- anharmonicity- 
mediated manipulation of THz polaritons by means of a bulk TA wave, is evaluated for 
CuCl as m 3 lU 2 - 0.26 meV (63 GHz) for J ac = lkW/cm 2 . In thallium halides (T1C1 and 
TIBr), quartic anharmonicity is dominant with positive values of V4. With the Leibfried- 
Ludwig approximation, one can evaluate the control parameter m^I^ from the available 



experimental data on the real part of the TO-phonon self-energy 18]. For T1C1 driven by a 
bulk TA wave we get m 4 / ac ~ 0.4 meV (97 GHz) for 7 ac = 100kW/cm 2 . Note that ra 3 and 
7714 scale to the sound velocity as ^ s 3 ^ 2 and v~ 3 , respectively, so that the use of surface AWs 
considerably reduces the operating acoustic intensity. 
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We examine the optical response of an acoustically-driven TO-phonon polariton in a one- 
dimensional geometry, when a semiconductor occupies the half space z > 0, and a normally 
incident light field of frequency uj induces a THz polariton propagating collinearly to the 
pumping AW (see Fig. 1). In this case, apart from reflectivity at the same frequency uj ) down- 
converted Bragg replicas at uj + tiuj^ (n = — 1, —2, ...) arise in the reflection spectrum, due to 
acoustically-induced backward scattering of the polariton. In order to calculate the multiple 
Bragg replicas, we develop an approach more advanced than that used in conventional 
acousto-optics. For the latter, the acousto-optical susceptibilities are so weak that usually 
only one Bragg replica n = +1 or —1 is seen. In contrast, the TO-phonon resonance 
mediates and considerably enhances the coupling between the optical and acoustic fields, 
so that generally one has to take into account the whole series of the Bragg replicas and 
multi-phonon transitions, thus treating the problem non-pert urbatively. 

The quasi-energy polariton spectrum uo = o;(p), calculated for real- valued p (quasiparticle 
solution) by solving Eqs. (E])-((3]) with = 0, is plotted in Fig. 2 (a) for CuCl parametrically 
driven by the TA wave of frequency z/ ac = 50 MHz and constant intensity 7 ac = 25kW/cm 2 . 
The spectrum, which arises from spatial and temporal modulation of the crystal lattice, 
can be interpreted in terms of the Brillouin zone picture, with acoustically-induced energy 
gaps A 7v oc I&J 2 due to the A-phonon resonant transitions within the polariton dispersion 
branches. The spectral positions of the transitions are indicated in Fig. 2 (a) by the vertical 
arrows. For the frequency scale used in Fig. 2(a), only the gaps Ajv=i and A N=2 in the 
upper polariton branch are clearly seen. 

From the air side, z < (see Fig. 1), the light field is given by 

E(z < 0, t) = e iqoz e~ iujt + rn6 -^ e -^+na; ac )t ^ ^ 

n 

where q n = (uj + nuj^ c )/c with — n max ^ n ^ n max (we proceed up to n max = 60) and r n 
stands for the amplitude of the outgoing Bragg replica n normalized to the unity amplitude 
of the incoming light wave. The electric field propagating in the crystal [z > 0) is 

E(z > 0, t) = AjEnje^^-^™^ . (5) 

Here, pj = Pj(cj) is the wavevector associated with the quasi-energy dispersion branch j 
(n max ^ j ^ n max ) of the acoustically-driven polariton, E n j are the corresponding normalized 
eigenvectors, and Aj are the eigenmode amplitudes. Both pj = Pj{uS) and E n j are the 
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FIG. 2: (Color online) (a) The dispersion uj — uj(p) of THz polaritons in CuCl driven by the 
bulk TA-wave of ^ ac = 50 MHz (black lines) and the bare polariton spectrum (red lines). The 
wavevector p is normalized to the acoustic wavevector k. The arrows highlight the 7V-TA-phonon 
resonant coupling between the polariton states (N = 1,2). (b)-(e) The calculated total Bragg 
reflectivity (black lines), R = J2 n \r n \ 2 , against the light frequency v — cj/(2tv) for ^ ac = 50 MHz, 
100 MHz, 200 MHz, and 1GHz, respectively. The bare polariton reflectivity is shown by the red 
lines. The Bragg signals labelled in (c) as n — — 1 and n — — 2 are mainly due to |r_i| 2 and |r_2| 2 , 
respectively. J ac = 25kW/cm 2 , htt TO = 20.28 meV, htt R = 14.53 meV, and ^7x0 = 0.2 meV. 

forced-harmonic solutions of Eqs. (ED-® for real- valued frequency uj. The exponential on 
the r.h.s. of Eq. ([5j) as well as the basic relationships pj +s (u)) = pj(uj — suj^ c ) + sk and 
E n j+ s (oj) = E n + S j(oj — scc; ac ), with integer 5, reflect the acoustic wavevector and frequency 
translational invariance of the quasi-energy spectrum. The Maxwellian boundary conditions 
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FIG. 3: (Color online) (a) The interaction length between the THz polariton and TA pumping 
wave in CuCl, Ant = Ant(^) 3 calculated for I ac = 25kW/cm 2 and z/ ac = 50 MHz, 100 MHz, and 
200 MHz. (b) The electric field profile \E{z > 0,cj)| 2 evaluated for = 100 MHz (A ac = 20.2 /mi). 
The color scale is logarithmic with black color corresponding to E — > 0. 



at z = together with Eqs. (Sj)-([5j) yield: 

^n,0 + r n = AjE n j , 



(6) 



The set of 2(2n max + 1) linear Eqs. ([6j) determines r n and Aj. 

In Figs. 2 (b)-(e) we compare the calculated total Bragg reflectivity, R = R(uj) = \r n \ 2 
(black solid lines), with the reflectivity of the acoustically unperturbed THz polariton, = 
R(°'(lj) (red solid lines). The sharp spikes in the Bragg spectrum of the acoustically-driven 
CuCl crystal are clearly seen for the one-TA-phonon and two-TA-phonon transitions, both 
for the upper (UP) and lower (LP) polariton branches [see Fig. 2 (a) against Fig. 2 (b)]. For 
is^c ^ 30 — 300 MHz, the backward scattered Bragg replica |r_i| 2 peaks at the energy of 
the one-phonon transition and is highly efficient, with |r_i| 2 /i? ~ 50% — 70%. The peak 
position and its strength are effectively tunable by changing the frequency and intensity 
of the AW [see Figs. 2 (b)-(d)]. This can be used for the frequency down- conversion by 
cj ac = 27rz/ ac of the optically-induced THz polariton and of the incident light field. Generally, 
the backward Bragg scattering signals |r n | 2 oc /ic' (n < 0) peak at the spectral position 
of the gaps A^^| n |, i.e., for the light frequency uu = uun which satisfies the resonant Bragg 
condition Pj=o{(jJn) — Nk/2 with N — 1, 2, .., |n|. 

The Bragg signal n = appears as the AW-induced change of the reflectivity at incident 
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frequency u ) |r (cc;)| 2 — R(°\oj). The strength of the n = signal sharply increases with 
decreasing detuning \uj — ^to| from the TO-phonon resonance, so that |ro| 2 — becomes 
dominant over |r_i| 2 . For z/ ac > 1GHz, when the one-acoustic-phonon transition within 
the LP branches occurs very close to f^ro? the AW-induced change of the reflectivity [see 
Fig. 2 (e)] is completely determined by the n = replica and has no z/ ac -down-converted 
frequency components. Thus this operating mode can be used for TO-phonon polariton 
deflectors and acoustically-controlled THz filters. 

The interaction length £- int = Ant(^) required for the formation of the Bragg signals and 
thus for AW control of the THz light field is plotted in Fig. 3 (a) for various i/ ac . The sharp 
troughs at uj = ujn=i and ujn=2 in the £- int = £- m t(uj) profile are due to the n = — 1 and —2 
Bragg replicas [see Fig. 3 (a)]. Figure 3 (b) shows the light field distribution associated with 
the optically-induced THz polariton in the acoustically-driven CuCl crystal. Apart from 
the broad black band [see Fig. 3 (b)], which corresponds to the Restrahlen band with rather 
weak penetration of the light field into the crystal, the narrow stripes of alternating color 
illustrate the formation of the Bragg replicas. The interaction length for the n — — 1 replica 
at its resonant frequency uj = ujn=i is given by 

Equation ([7]), which is valid for |^to — ^jv=i| ^> 7to and £- mt k 3> 1, shows the resonant 
decrease of £- int oc 1/ vXc with decreasing frequency detuning from the TO-phonon resonance. 
In this case the interaction between the light field and pumping AW is mediated by the TO- 
phonon resonance, giving rise to £- int of only a few tens of acoustic wavelength A ac (see Fig. 3). 
This is in sharp contrasts with conventional acousto-optics where £- int ~ 10 3 — 10 4 A ac for the 
same operating I ac . 

Within the used nonperturbative approach, each Bragg replica n integrates all n + s — s 
TA-phonon transitions with s ^ n max : With increasing 7 ac the bare n-phonon transitions 
become dressed by higher-order processes when n + s phonons are emitted and s phonons 
absorbed. For the n = — 1 replica shown in Fig. 2(b), e.g., the multiphonon transitions 
— 1 + 1 — 1, —1 + 2 — 2, etc. account for about 90 % of |r_i| 2 . For v ac > 1 GHz, the dominant 
contribution to |ro| 2 — stems from —s + s multi- TA-phonon transitions with s > 1. 

The calculated room temperature reflectivity R = R(uj) of a T1C1 crystal driven by the TA 
wave of frequency u ac = 25 MHz and 125 MHz is plotted in Figs. 4 (a) and (b), respectively. 
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FIG. 4: (Color online) The total Bragg reflectivity R = R(ou) of a T1C1 crystal driven by the TA 
wave of frequency u ac = 25 MHz (a) and 125 MHz (b), I ac = 100kW/cm 2 (black dotted lines) 
and 200kW/cm 2 (black solid lines). The acoustically-unperturbed THz spectrum = R(°\oj) 
is shown by the red lines. Insert: The acoustically-induced Stark shift of the TO-phonon line. 
Ml T o = 7.81 meV, hQ R = 17.97 meV, and hyro = 0.92 meV. 

In this case, in Eqs. (pQ) and (J3J) we put m 3 = 0, and only even-order TA-phonon-assisted 
transitions occur. The Bragg signals n = —2, due to |r_ 2 | 2 oc 7f c , are indicated in Fig. 4 (a) 
by arrows, for the transitions within the LP and UP branches, respectively. Similarly to the 
previous case (CuCl), for THz light frequency uu close to Qto the AW-induced change of i?, 
AR — R — R(°) ~ 7 2 c [see Fig. 4 (b)], is mainly due to the n = Bragg replica. The quartic 
nonlinearity leads to the Stark blue shift by 2m^l^ of the TO-phonon frequency, according 
to Eqs. (PQ) and (Ej), as is clearly seen in Figs. 4 (a) and (b). The Stark shift ~ 0.1 - 0.2 THz 
has a rather sharp contrast on the blue side of the THz reflectivity [see inset in Fig. 4 (b)]. 

The acoustically-induced modulation of infrared polaritons has to be particularly strong 
br ferroelectric soft TO-phonons (e.g., in LiTa03 and LiNb03 [6] and bismuth titanate 
91). In this case, a multi-well local potential for the displacive ferroelectric mode has a 
considerable low-wavevector component and therefore yields large values of V3 and V4. Far- 
infrared optical phonons (2— lOmeV) in zirconium tungstate (ZrW^Os) indicate anomalously 



high anharmonicity 



, The normal modes associated with soft TO-phonons in this 
negative thermal expansion compound are a mixture of librational and translational motion. 
The latter strongly couples with acoustic phonons giving rise, as we foresee, to manipulation 
of the THz polaritons by using utrasound waves of modest I ac . 
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